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1.  Introduction 

Three  decades  of  incoherent  scatter  radar  measurements  have  greatly  advanced  ionospheric 
physics  by  providing  accurate  electron  and  ion  densities,  electron  and  ion  temperatures,  and  ion 
velocities  as  a  function  of  time  and  altitude.  Reviews  on  incoherent  scatter  radar  measurements 
and  techniques  were  made  by  Burnside  et  al.  (1991),  Suvanto  (1990),  Mathews  (1986),  Sulzer 
(1986),  Evans  (1978),  and  Farley  (1971).  By  now  there  are  long  term  datasets  in  most  of  the 
incoherent  radar  observatories,  harboring  what  can  be  termed  as  ionospheric  climatology.  This 
report  presents  frequency  histograms  covering  ten  years  of  incoherent  scatter  radar  measurements 
at  the  Arecibo  Observatory  (18°  20’  N,  66°  45’  W,  50°  magnetic  dip)  in  Puerto  Rico.  The 
Arecibo  measurements  were  made  during  the  regularly  scheduled  World  Days  using  the  multi¬ 
frequency  technique  of  Sulzer  (1986)  between  October  1985  and  May  1995  [Figure  1].  Over 
140,000  Arecibo  profiles  were  recently  reanalyzed  at  the  Observatory  in  order  to  benefit  from 
improvements  in  numerical  fitting  of  the  measured  incoherent  scatter  spectra.  These 
improvements  led  to  resolving  0%  H^,  and  He'^  ions  simultaneously  in  the  measured  spectra 
[Gonzalez  and  Sulzer,  1996],  and  retrieving  more  accurate  topside  ion  and  electron  temperatures 
[Sulzer  and  Gonzalez,  1996]. 

The  standard  World  Day  mode  at  Arecibo  records  incoherent  scatter  spectra  for  altitudes  between 
144  and  627  km,  in  38-km  increments.  Before  1985,  the  standard  World  Day  program 
transmitted  300  ps  pulses  at  430.0  MHz.  Since  October,  1985,  the  program  uses  the  multi¬ 
frequency  technique  of  Sulzer  (1986),  phase  modulating  a  308-ps  pulse  and  recording  seven 
independent  spectra  per  range  gate.  These  are  alternated  with  a  13-baud  (4  ps  per  baud)  Barker 
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coded  pulse.  In  addition,  15  range  gates  per  IPP  (inter-pul.se  period)  are  recorded  after  September 
1986,  instead  of  14,  thanks  to  a  faster  array  processor. 

World  Day  observations  between  1985  and  1990  are  discussed  by  Burnside  et  al.  (1991).  Their 
14-day  period  averages  show  solar  cycle  control  of  the  height  of  the  F2  layer  peak  height  (zmax) 
and  density  (Nmax).  Other  studies  summarize  Arecibo  ExB  drifts  [Fejer,  1996],  electric  fields 
[Ganguly  et  al.,  1987],  and  meridional  wind  data  [Harper,  1973].  This  study  presents  a  statistical 
summary  or  climatology  of  the  Arecibo  World  Days  in  the  NCAR-CEDAR  database  up  to  1995. 
Specifically,  we  present  frequency  distributions  of  important  observables  such  as  peak  density, 
peak  height,  electron  temperature  (Te),  the  electron  and  ion  temperature  ratio  (T/Ti),  line-of-sight 
ion  velocities,  and  and  He'^  fractions.  No  published  report  we  know  shows  long-term 
statistical  frequencies  of  Arecibo  measurements.  The  importance  of  the  distributions  is  in 
quantifying  the  physically  realizable  values  of  the  ionospheric  state  variables  measured  during 
World  Days.  Such  information  is  also  valuable  in  evaluating  numerical  or  empirical  models, 
geophysical  events,  and  in  planning  and  designing  future  sensor  and  experiments.  Furthermore, 
the  Arecibo  electron  density  profiles  are  fitted  using  a  four-parameter  model,  and  the  parameter 
distributions  are  presented. 


2.  The  “Linear-H”  Model 


A  model  of  the  electron  density  profile  is  required  to  extract  F2  jjeak  height  and  density  from  the 
measured  altitude  profiles  of  electron  density.  This  is  because  the  relatively  coarse  altitude 
resolution  (~38  km)  may  miss  the  actual  peak  height.  In  addition,  noise  can  obscure  the  precise 
values  of  both  the  peak  height  and  density.  In  the  World  Day  program,  absolute  power 
calibration  is  provided  by  scaling  the  density  profiles  with  foF2  measurements  from  the 
ionosonde  on-site  or  the  Ramey  ionosonde  some  40-km  to  the  northwest.  However,  the 
ionosonde  calibration  is  less  frequent  than  the  measurement  interval  of  a  few  minutes  per  profile. 
The  simplest  way  to  extract  the  parameters  at  the  peak  is  just  to  define  the  altitude  with  the 
largest  electron  density  as  zmax,  with  Nmax  as  the  density  at  that  height.  The  obvious  drawback 
is  that  the  actual  peak  may  be  located  above  or  below  the  range  gate  of  maximum  electron 
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density.  A  popular  approach  is  to  fit  with  a  function,  such  as  a  parabola,  the  three  consecutive 
range  gates  centered  on  the  largest  electron  density.  The  zmax  and  Nmax  reported  in  the  NCAR- 
CEDAR  database  are  obtained  by  fitting  the  entire  profile  with  cubic  splines,  an  interpolation  that 
goes  through  each  data  point.  These  methods  are  flawed  in  that  they  do  not  account  for  the  ever¬ 
present  noise,  giving  x^=0.  However,  there  is  generally  good  agreement  between  the  cubic- 
splines  fitting  and  the  Barker-coded  profiles.  The  Barker-coded  profiles  have  600-m  height 
resolution. 


Our  approach  is  to  fit  the  measured  electron  density  profile  with  a  four-parameter  Chapman-like 
function  [Picone  et  al.,  1997]  or  “linear-H”  model,  excluding  range  gates  below  200  km.  The 
Chapman  layer  representation  [Chapman,  1931],  and  an  altitude  varying  scale  height  are 
discussed  by  Rishbeth  and  Garriott  [1969].  Other  Chapman-like  layers  are  discussed  by  Fox 
[1994].  We  denote  the  model  peak  height  and  density  as  hmax  and  nmax,  respectively,  to 
distinguish  them  from  the  corresponding  CEDAR  database  parameters  zmax  and  Nmax, 
respectively.  The  function  is  similar  to  the  three-parameter  Chapman  formula  [Burnside  et  al, 
1983;  Anderson  et  al.,  1985]  but  adds  a  fourth  parameter,  which  defines  the  height-varying 
pseudo-scale  height  linear  term,  hence  the  name  linear-H.  The  linear-H  function  is  given  by 


Ne(z) 


nmax  exp 


(1) 


where  z(h)  is  the  reduced  height  given  by 

z(h)  =  (h  -  hmax)  /  Hm(h)  ( 1  a) 

and  Hm(h)  is  the  linear  scale  height  parameter,  given  by 

Hm  =  Ho  +  Hi  (h  -  hmax)  for  h  >  hmax  and  Hm  =  Ho  for  h  <  hmax.  (lb) 

The  four  parameter  vector  for  the  model  is  thus  [nmax,  hmax.  Ho,  H|].  Numerical  fits  are  made 
using  Levenberg-Marquardt  minimization  of  [Picone  et  al,  1997],  which  accounts  for 
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measurement  noise.  Two  advantages  of  the  linear-H  function  are  simplicity  —only  four 
parameters  are  used  to  describe  the  profile—  and  that  the  parameters  are  geophysically 
meaningful.  The  physical  basis  of  the  model  is  that  the  function  is  a  solution  to  the  ion  continuity 
equation.  Below  we  provide  statistical  measures  of  the  quality  of  the  resulting  fits. 


3.  Results 

3.1  Dataset 

Table  1  summarizes  the  statistical  properties  of  the  linear-H  fits  to  the  Arecibo  World  Days 
between  1985  and  1995.  The  distribution  of  the  daily  F10.7  or  10.7-cm  flux  values 
corresponding  to  the  Arecibo  World  Days  considered  here  is  shown  in  Figure  1.  High  solar 
activity  is  centered  in  the  middle  of  the  observing  period  during  1989-1992,  with  F10.7  index 
reaching  275.  Low  solar  activity  is  found  at  the  beginning  and  ending  periods,  with  F10.7  values 
as  low  as  75  during  1985-87  and  1994-95.  Since  the  median  F10.7  is  133,  the  dataset  contains 
more  World  Days  in  solar  maximum  than  solar  minimum.  This  small  bias  results  from  a  data  gap 
in  the  second  half  of  1991,  as  shown  in  Figure  2.  Monthly  coverage  is  irregular,  as  shown  in 
Figure  3:  March,  May,  and  July  are  barely  represented  in  the  dataset  with  frequencies  of  less  than 
0.05  or  5%  each.  Most  of  the  measurements  were  taken  during  the  months  of  January,  April, 

July,  and  during  Fall.  Local  time  coverage  is,  on  the  other  hand,  highly  uniform,  as  shown  in 
Figure  4.  Antenna  pointing,  shown  in  Figure  5,  is  very  uniform  with  a  beam  elevation  of  15°  off- 
zenith.  As  is  typical  of  the  World  Day  program,  the  beam  is  swung  azimuthally,  recording  three 
profiles  every  four  minutes  during  the  1 6  minutes  needed  for  one  360°  rotation. 

The  distribution  of  peak  densities  in  Figure  6  shows  large  positive  skew,  with  the  most  frequent 
bin,  denoted  by  “mode”,  at  0.2  x  lO'^  m'^  (Table  1).  Two  histograms  are  shown:  one  for  the 
nmax  obtained  with  linear-H,  another  for  the  Nmax  values  directly  from  the  NCAR-CEDAR 
database  (obtained  with  cubic  spline  interpolation).  The  linear-H  and  splines  histograms  are  very 
similar,  though  the  CEDAR  Nmax  has  mean  and  median  which  are  20%  and  30%  smaller  than 
the  corresponding  linear-H  values  (Table  1).  The  reason  for  this  is  that  the  linear-H  model  has 
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distribution  has  slightly  lower  nmax  values  below  about  0.4  x  lO'^  but  a  higher  frequency  of 
larger  peak  densities  than  CEDAR  between  0.4  x  10*^  m'^  and  2  x  10*^  m'^.  Only  23  out  of  9801 
fits  are  larger  than  3x10  m  ‘ .  The  distribution  includes  both  day  and  night  measurements  with 
the  largest  peak  densities  occurring  in  the  afternoon  [Burnside  et  al.,  1991].  Using  data  from  half 
of  a  solar  cycle,  Burnside  et  al.  (1991)  found  Nmax  to  increase  with  increasing  F10.7  values.  Our 
data  shows  a  similar  trend. 

The  peak  height  distributions  in  Figure  7,  on  the  other  hand,  show  a  broad  and  symmetric  linear- 
H  hmax  distribution  (solid  line)  with  mode  at  315  km,  while  the  corresponding  CEDAR  (splines) 
distribution  fluctuates  considerably.  The  two  distributions  have  similar  statistics  (Table  1).  The 
mean  hmax  is  316  ±  51  km  (“mean  ±  standard  deviation”)  for  the  linear-H  model;  the  CEDAR 
mean  is  320  ±  49  km  (Table  1).  The  median  in  the  linear-H  distribution  is  3 17  km  whereas  the 
median  in  the  CEDAR  hmax  is  316  km.  The  CEDAR  distribution  has  two  large  frequencies  at 
about  310  km  and  365  km.  The  linear-H  distribution  can  be  fitted  with  a  Gaussian  function  of 
amplitude  0.02,  mean  of  314.3  km,  and  standard  deviation  of  52.3  km.  These  parameters  agree 
with  those  in  Table  1  to  better  than  3  km,  suggesting  that  the  linear-H  hmax  is  normally 
distributed  unlike  nmax  (Fig.  6).  The  Gaussian  is  most  accurate  for  hmax  values  above  about  370 
km  while  it  overestimates  the  very  low  peak  heights. 

The  highest  F2  layer  at  Arecibo  occurs  prior  to  midnight,  reaching  380  km  [Burnside  and  Tepley, 
1989].  The  phase  (i.e.,  local  time)  of  highest  hmax  is  also  prior  to  midnight  with  the  linear-H  fits 
to  the  data.  The  lowest  layer  height  occurs  in  the  morning,  reaching  an  average  of  260  km 
[Burnside  et  al.,  1991].  The  extremes  in  the  diurnal  average  of  Burnside  et  al.  fall  within  the 
standard  deviation  of  the  distribution.  Our  distribution  is  consistent  with  the  scatter  in  the 
diurnal  variation  reported  by  Burnside  et  al.,  which  can  be  as  large  as  100  km.  They  also  found 
solar  cycle  variability  to  be  larger  than  the  seasonal  variability.  Our  preliminary  evaluation  of  the 
1985-1995  fits  gives  results  similar  to  those  of  Burnside  et  al.  Our  database  has  the  advantage  of 
greater  sampling,  particularly  during  high  F10.7  values.  Thus,  the  distributions  in  Figure  7  are 
dominated  by  diurnal,  solar  cycle  and  seasonal  variations,  in  descending  order  of  importance. 
Overall,  the  Arecibo  F2  peak  is  rarely  found  above  400  km  or  below  230  km  though  peak  heights 
below  200  km  are  less  accurately  determined.  This  is  due  to  the  large  bottomside  gradient  in  the 
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electron  density  caused  by  the  higher  loss  rates,  also  known  as  the  range-smearing  problem  that 
occurs  when  the  pulse  length  is  longer  than  the  characteristic  scale  length  of  the  medium. 

3.2.  Model  and  Data  Uncertainties 

A  model  fit  should  provide,  in  addition  to  the  fit  parameters  and  uncertainties  (“errors”),  a 
statistical  measure  of  the  goodness  of  fit.  The  distribution  of  the  reduced-^^  values  obtained  in 
fitting  the  profiles  with  the  linear-H  model  is  shown  in  Figure  8;  the  calculation  includes  both 
measurement  and  model  error  covariance  calculations  as  described  by  Picone  et  al.  [1997]. 

Figure  8  shows  the  distribution  of  logioof  reduced  one  value  per  fitted  profile.  The  fits  to  the 
chosen  Ne(h)  profiles  are  obtained  by  minimizing  with  a  modified  Levenberg-Marquardt 
algorithm,  as  discussed  by  Picone  et  al.  [1997].  Each  profile  fit  is  terminated  when  changes  by 
less  than  1  part  in  10'^  or  upon  reaching  200  iterations,  whichever  comes  first.  The  distribution 
peaks  at  logiox^  =  0  or  =  1,  indicating  a  good  overall  fit.  A  value  of  x^  =  1  means  that  the 
residuals  or  deviations  are  equal  to  the  random  fluctuations  expected  from  the  measurements. 

The  overall  distribution  has  a  mean  of  x^=  1.15  ±  6.09  (logioX^  =  0.06  ±  0.785).  This  is  the 
result  of  the  model  parameters  chosen  by  the  fitting  algorithm,  excluding  fits  with  Ix^l  <  lO  '*. 
Values  of  x^  greater  than  10  and  smaller  than  0.1  account  for  less  than  16%  of  all  the  fits  (Table 
1);  thus,  16%  represents  the  nominal  fraction  of  “outlier”  fits.  Therefore,  most  of  the  fits  produce 
adequate  values  of  y}.  Too  small  x^  fits  may  indicate  large  or  overestimated  experimental 
uncertainties,  while  large  y  fits  may  indicate  underestimated  data  uncertainties  or  a  bad 
functional  fit. 

Another  measure  of  the  accuracy  of  the  model  fits  is  the  distribution  of  the  uncertainties  in  the 
data  and  fit  parameters.  By  uncertainty,  we  mean  a  variety  of  statistical  measures  of  error,  such 
as  the  residual  {x  -  X)  between  the  model  (linear-H)  parameter  jc  and  the  corresponding  data 
parameter  X.  Also,  AX  denotes  the  uncertainty  in  the  CEDAR  data  while  Ax  denotes  the 
uncertainty  in  the  model  (obtained  from  the  diagonal  element  of  the  covariance  matrix  of  the  fit), 
normalized  by  either  x  or  X.  The  quantity  a  denotes  the  standard  deviation,  specifically  0  = 
V(variance). 
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Figure  9  shows  the  relative  model  fit  uncertainty  in  nmax,  specifically,  the  fractional  percentage 
uncertainty  in  the  model  peak  density,  Anmax,  given  by 

Anmax  =  100.  a(nmax)  /  nmax  (2) 

where  a(nmax)  is  obtained  from  the  linear-H  fit  covariance  matrix.  The  covariance  matrix 
equation  is  given  by  equation  (13)  of  Picone  et  al.  [1997],  and  the  variance  is  estimated  using  a 
formula  similar  to  Eqn.  (15)  from  that  reference.  The  overall  distribution  looks  like  the  sum  of  a 
Gaussian  and  a  smaller  skewed  distribution  (Fig.  9).  The  skewed  distribution  is  not  negligible 
since  nearly  28%  of  the  distribution  has  values  above  4%  though  the  low  frequencies  (below 
0.002)  indicate  considerable  spread.  However,  85%  of  the  Anmax  distribution  is  within  10% 
(Table  1).  The  “outlier”  population  (remaining  15%)  is  of  the  same  size  as  the  outliers  in  the 
nmax  and  hmax  distributions  (Table  1).  The  mode  (most  frequent  value)  in  the  distribution  of  the 
uncertainty  in  nmax  (Equation  2)  is  nearly  2%  (Fig.  9)  with  a  frequency  just  below  2%  also,  using 
200  bins.  There  are  no  negative  values  since  a  (nmax)  is  positive  definite.  The  median,  which  is 
less  sensitive  to  outliers,  is  2.4%.  If  Anmax  >  100%  are  excluded  (13%  of  the  distribution)  the 
mean  becomes  3%  ±  5%.  The  precision  in  the  model  peak  density  is  comparable  to  the  typical 
uncertainty  of  the  Arecibo  incoherent  scatter  radar  measurements  of  the  peak. 

The  distribution  of  the  uncertainties  in  the  model  hmax,  Ahmax,  shown  in  Figure  10,  is  analogous 
to  Anmax.  However,  Ahmax  is  about  three  times  narrower  distributed  than  Anmax,  and  its  mode 
is  twice  as  frequent  (Table  1).  The  median  in  Ahmax  is  0.6%,  or  about  4  times  smaller  than  the 
median  in  Anmax.  Both  the  Ahmax  and  Anmax  distributions  have  similar  populations  with 
values  of  10%  or  less  (Table  1).  However,  Ahmax  is  less  skewed  than  Anmax.  Neglecting 
Ahmax  >  100%  (12%  of  the  population),  the  mean  becomes  1.4%  ±  5.5%.  Therefore,  the  linear-H 
uncertainties  are  lower  for  hmax  than  for  nmax.  In  other  words,  the  data  provides  a  more  precise 
determination  of  hmax  than  nmax.  The  model  uncertainty  in  hmax  is  uniformly  distributed  in 
time  while  Anmax  peaks  after  dawn. 
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The  differences  between  the  linear-H  and  the  splines  models  can  be  described  by  comparing  the 
distribution  of  the  residuals  between  the  respective  values  of  nmax  and  hmax.  Figures  1 1  and  12 
show  the  distributions  of  the  percentage  residuals,  D,  in  peak  density  and  height,  respectively,  as 
given  by 

D(x-X)=  100.  (x-X)/X  (3) 

where  x  =  linear-H  fit  of  hmax  or  nmax,  and  X  =  CEDAR  (cubic  splines)  zmax  or  Nmax, 
normalized  by  the  corresponding  CEDAR  values.  Zero-residual  occurs  when  the  linear-H  and 
splines  estimates  are  equal.  The  asymmetry  and  negative  skew  in  the  peak  density  residuals  (Fig. 
11)  show  a  definitive  bias  between  the  two  estimates.  In  particular,  the  linear-H  peak  density  is 
lower  than  the  splines  estimate  in  58%  of  all  cases,  a  result  also  implied  in  Fig.  6.  The  bias  or 
mean  residual  is  -0.9  ±  4.5%;  the  median  is  also  small  and  negative  (Table  1).  On  the  other 
hand,  the  two  methods  usually  yield  similar  estimates.  Values  of  D(nmax-Nmax)  fall  outside 
±5%  in  less  than  10%  of  all  cases,  and  outside  ±10%  in  less  than  3%  of  all  cases  (Table  1).  In 
other  words,  the  great  majority  of  linear-H  and  the  CEDAR  database  estimates  are  within  10%  of 
each  other.  The  “3a”  level  of  confidence  of  a  normal  distribution  is  at  the  residual  level  of  ±10% 
since  these  account  for  nearly  97%  of  all  comparisons. 

The  parameter  D(hmax-zmax)  is  the  residual  between  the  model  and  measured  hmax,  and  is  more 
normally  distributed  (Fig.  12)  than  D(nmax-Nmax)  (Fig.  11).  There  is  a  negative  bias  of-1%, 
which  implies  that  the  linear-H  peak  height  is  three  km  lower  than  the  splines  on  average  (Table 
1).  However,  peak  height  estimates  by  the  two  methods  are  randomly  distributed  about  the  -1- 
km  mean.  Over  90%  of  the  values  of  D(hmax-zmax)  are  5%  or  smaller.  Therefore,  the  data  and 
the  linear-H  fits  agree  to  within  5%  in  well  over  90%  of  the  cases,  well  within  the  “2a”  or  90% 
level  of  confidence  for  a  normal  distribution.  The  linear-H  estimate  of  hmax  is  somewhat  more 
precise  and  accurate  than  its  estimate  of  Nmax  (Figs.  1 1-12).  There  is  also  a  negative  skew  in  the 
distribution,  where,  for  example,  the  -5%  residual  is  over  six  times  more  frequent  than  the  ±5% 
residual.  The  hmax  (linear-H)  is,  therefore,  more  often  smaller  than  zmax  (CEDAR),  a  finding 
supported  by  Fig.  7.  This  difference  peaks  at  mid-moming,  so  the  linear-H  value  is  smallest  at 
that  time  (not  shown).  The  model  value  is,  on  the  other  hand,  largest  at  midnight. 
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The  linear-H  model  is  actually  a  representation  of  the  electron  density  profile,  not  just  the  peak 
height  and  density.  A  measure  of  the  overall  accuracy  of  the  linear-H  fit  to  the  measured  electron 
density  profile  is  the  average  of  the  altitudinal  sum  of  the  absolute  deviations  between  the  model 
and  measured  electron  densities,  normalized  by  the  measured  density.  This  average  residual,  B,  is 
given  by 

B(ne)  =  100  (1/M)  Z  KneO)  -  Ne(j))l  /  NeO)  (4) 

where  ne(j)  is  the  linear-H  electron  density  at  the  jth  range  gate,  Ne(j)  is  the  corresponding 
measured  value,  and  M  is  the  number  of  range  gates  (altitudes)  in  each  profile.  The  sum  in  (4)  is 
over  all  range  gates  fitted  (i.e.,  at  or  above  200  km),  and  is  always  positive.  The  parameter  B  is 
the  average  magnitude  of  the  residual  of  the  model  and  measured  electron  densities  in  each  fitted 
electron  density  profile.  The  closer  B  is  to  zero,  the  more  similar  are  the  model  and  measured 
Ne(h)  profiles.  The  quantity  B  measures  how  well  a  profile  was  fitted  while  the  distribution 
measures  how  many  profiles  were  fit  with  a  particular  B  value.  The  distribution  of  B,  shown  in 
Figure  13,  has  a  mode  between  1%  and  2%,  and  a  plateau  centered  at  5%.  The  mean  and 
standard  deviation  in  the  distribution  of  5  are  5  ±  23%,  and  the  median  is  3.6%  (Table  1).  The 
relatively  large  standard  deviation  is  caused  by  the  positive  skew  in  the  distribution,  and  a  few 
large  outliers.  However,  the  cumulative  frequency  for  B  >  20%  is  less  than  0.015  (1.5%).  That 
is,  less  than  2%  of  all  fitted  profiles  have  an  average  residual  of  more  than  20%  (Table  1).  If  the 
distribution  is  truncated  to  B  <  100%,  the  mean  is  unchanged  but  the  standard  deviation  decreases 
to  ±  5%.  Inspection  of  Fig.  13  shows  that  over  69%  of  the  values  of  B  are  smaller  than  7%. 

Since  94%  of  all  fits  have  values  of  B  <  10%,  the  linear-H  model  fits  the  measured  profiles  with 
an  accuracy  of  better  than  10%  at  the  “2a”  or  90%  level  of  confidence  for  a  normal  distribution. 
The  best  fits,  as  measured  with  B,  are  achieved  at  night  whereas  the  values  of  B  near  6%  occur  in 
the  morning  (not  shown).  As  expected,  the  linear-H  model  is  most  accurate  at  night  when  there  is 
no  E  or  Fi  layer. 
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There  is  an  analogous  parameter  to  B,  but  one  based  solely  on  the  uncertainties  of  the  measured 
electron  density  profile  listed  in  the  CEDAR  database.  This  is  the  altitude-averaged  percentage  of 
the  experimental  uncertainty  in  Nc(h)  normalized  by  the  measured  electron  density,  D(Ne),  given 
by 


D(Ne)  =  100  (1/M)  Z  a(Ne(j))  /  KQ)  (5) 

where  c(Ne(j))  is  the  uncertainty  in  the  measured  Ne(j)  at  the  jth  range  gate,  summed  over  all  the 
range  gates;  the  other  quantities  are  as  defined  in  equation  (4).  Thus,  D(Ne)  measures  the  average 
experimental  uncertainty  of  all  profiles.  This  quantity  is  readily  compared  to  the  average  fit 
residual  B  (Fig.  13).  Similar  distributions  of  B  and  D  would  indicate  similar  levels  of  precision 
in  the  linear-H  and  measured  electron  density  profiles.  Such  comparison  (Table  1)  indicates  that 
the  accuracy  of  the  linear-H  model  is  similar  to  the  average  experimental  uncertainty  since  the 
mean,  median,  and  standard  deviation  of  B  and  D  are  nearly  identical  (Table  1).  That  is,  the 
model  fit  is  as  accurate  as  the  data  is  precise.  The  distributions  are  different  in  that  D  is  nearly 
Gaussian  whereas  B  is  skewed. 

The  two  distributions  have  identical  cumulative  frequencies;  over  98%  of  the  distributions  are 
20%  or  less  in  value.  Less  than  10%  of  the  D  distribution  has  values  larger  than  10%,  a  figure 
comparable  to  that  of  the  outliers  in  the  model  fits.  The  D  distribution  has,  however,  more  and 
larger  outliers  than  the  B  distribution.  If,  on  the  other  hand,  profiles  with  D  >  100%  are  ignored, 
the  mean  becomes  5.6  ±  4.3%  which  is  just  about  equal  to  the  median  and  the  corresponding 
mean  of  the  B  <  100%  distribution  (5  ±  5%).  Therefore,  the  uncertainties  in  the  model  and 
measured  profiles  are  similar  though  their  distributions  have  different  shapes. 

The  distribution  of  the  difference  in  the  previous  two  measures,  B  -  D,  is  symmetric  and 
Gaussian-like,  as  seen  in  Figure  15.  This  can  be  interpreted  as  the  distribution  of  the  average 
difference  in  the  uncertainty  between  the  linear-H  and  measured  electron  density  profiles.  The 
difference  is  positive  when  data  errors  dominate,  and  negative  when  linear-H  errors  dominate. 
Values  greater  than  ±  100%  have  a  cumulative  frequency  of  less  than  0.0004.  The  mean  of  the 
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distribution  becomes  -1  ±  6%  if  only  differences  within  ±100%  are  counted,  with  median  of  - 
1.4%  (Table  1).  Therefore,  the  model  and  data  uncertainties  are  within  a  few  percent  of  each 
other,  with  slightly  larger  overall  data  errors.  Moreover,  98%  of  the  distribution  is  within  ±  20%, 
94%  is  within  ±  10%,  and  71%  is  within  ±  5%.  Less  than  5%  of  the  truncated  distribution  falls 
outside  of  ±  10%. 

Each  electron  density  profile  has  a  vector  of  absolute  residuals,  F,  normalized  by  the  data,  given 
by 


F(h)  =  100  l(ne(h)  -  Ne(h))l  /  Ne(h)  (6) 

where  n^fh)  is  the  linear-H  fit  to  the  measured  electron  density  at  height  h,  Ne(h).  The  median  of 
each  array  of  F  contributes  one  value  to  the  distribution  of  Figure  16.  The  median  of  F  may  be 
regarded  as  the  median  of  the  terms  in  the  sum  of  B  (Eqn.  4).  That  is,  B  =  ZF,  where  the  sum  is 
over  the  range  gates  of  a  given  profile.  The  most  frequent  value  of  F  is  smaller  than  1%,  and  the 
distribution  is  narrow  (Table  1).  Few  profiles  have  median  values  of  10%  or  higher;  their 
cumulative  frequency  is  less  than  0.05.  Over  25%  of  the  profiles  in  Fig.  16  have  a  median  F 
between  0%  and  1%.  The  distribution  of  medians  of  F  is  consistent  with  the  magnitudes  and 
mean  values  of  the  distributions  in  D  and  {B  -  D).  An  appreciable  number  of  large  median 
values,  not  seen  in  Figure  16,  would  denote  asymmetrical  or  strongly  skewed  distributions. 

The  minimum  uncertainty  of  each  measured  profile  as  a  fraction  of  the  linear-H  density,  C,  is 
given  by 

C  =  MIN  [a(Ne)  /  tie]  (7) 

where  o(Ne)  is  the  array  of  uncertainties  in  the  measured  electron  density  profile  (CEDAR 
database),  and  ne  is  the  linear-H  electron  density  profile  fit.  The  distribution  of  C,  shown  in 
Figure  17,  is  the  combination  of  a  near-zero  spike,  a  symmetrical  Gaussian-like  function,  and  a 
small  positive  skew.  The  spike  is  a  consequence  of  the  relatively  large  distribution  of  errors 
smaller  than  0.5%.  There  is  a  secondary  broad  peak  at  about  2%,  with  a  cumulative  frequency  of 
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0.05.  The  mean  of  the  distribution  is  2.2  ±  1.5%,  and  the  median  is  2%.  Moreover,  99.7%  of  C 
are  within  10%,  and  96.6%  is  within  5%.  Therefore,  the  most  frequent  minimum  experimental 
uncertainty  is  several  orders  of  magnitude  smaller  than  the  linear-H  density.  The  C  distribution  is 
narrower  than  that  of  the  average  profile  uncertainty  (Fig.  13).  Nearly  half  of  the  measured 
profiles  have  values  of  C  smaller  than  2%. 

The  statistical  measures  show  that  the  linear-H  model  is  accurate  to  within  5%  level  in  well  over 
67%  of  the  fitted  profiles,  and  accurate  to  within  10%  in  at  least  90%  of  the  profiles.  An 
examination  of  the  remaining  model  parameters  in  (1)  is  in  order.  Whereas  the  model  parameters 
nmax  and  hmax  define  the  magnitude  and  location  of  the  peak  electron  density.  Ho  and  Hi  define 
the  height  variation  or  profile  shape.  These  parameters  are  deemed  useful  in  and  of  themselves 
given  the  demonstrated  accuracy  of  the  linear-H  model  fits. 

3.3.  Model  Profile  Shape  Parameters  Ho  and  H| 

The  parameters  Ho  and  Hi  are  related  theoretically  to  neutral  temperature  and  its  altitude 
variation,  respectively,  though,  strictly  speaking,  they  are  numerical  parameters.  The  sum  Ho  + 

Hi  represents,  in  theory,  the  atomic  oxygen  scale  height  H  =  kT(h)  /  Mg,  where  M  is  the  mass  of 
atomic  oxygen,  T(h)  is  the  neutral  temperature  at  height  h,  and  g  is  the  acceleration  of  gravity. 
More  precisely,  Ho  +  Hi  is  an  effective  or  pseudo-scale  height  usually  taking  values  consistent 
with  an  atomic  oxygen  scale  height.  The  parameter  Ho  is  associated  with  the  scale  height  near 
the  F2  peak  while  Hi  is  the  average  height  gradient  above  the  peak  (Eqn.  1).  As  fit  parameters. 

Ho  and  Hi  can  have  whatever  values  are  assigned  to  them  in  order  to  minimize  while  matching 
the  shape  of  the  measured  profile.  The  parameter  Ho  is  on  the  order  of  50  km,  as  seen  in  Figure 
18.  This  value  corresponds  to  the  atomic  oxygen  scale  height  of  a  typical  thermospheric  neutral 
temperature  of  1000  K.  The  distribution  of  values  of  Ho  in  Figure  18  has  a  broad  peak  at  about 
50  km,  and  a  spike  at  20  km.  About  85%  of  the  Ho  values  forming  the  spike  at  20  km  are  for 
F10.7  values  of  100  or  less.  In  other  words,  the  spike  is  largely  due  to  narrow  profiles  during 
solar  minimum.  There  are  values  of  Ho  as  high  as  120  km  and  as  low  as  15  km,  for  a  six-fold 
spread.  The  mean  in  Ho  =  51  ±  16  km,  and  the  median  is  50  km  (Table  1).  The  extreme  values 
are  probably  less  representative  of  an  actual  scale  height. 
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The  spread  of  values  in  Fig.  18  represents  the  variability  of  electron  density  profile  shapes  since 
Ho  controls  the  altitudinal  change  of  the  linear-H  profile,  particularly  below  and  just  above  hmax. 
The  distribution  of  Ho  can  be  regarded  as  a  proxy  for  the  distribution  in  the  neutral  temperature, 

T,  since  the  scale  height  is  directly  proportional  to  T.  Thermospheric  temperature  increases  with 
increasing  F10.7  by  over  50%  [Hedin,  1987].  Since  the  distribution  covers  both  large  and  small 
FI 0.7  values  (Fig.  1),  some  of  the  larger  values  of  Ho  are  due  to  the  increased  solar  heating  of  the 
upper  atmosphere  during  solar  maximum  conditions.  The  average  time  variation  in  Ho  (not 
shown)  peaks  in  the  late  afternoon,  consistent  with  the  measured  variation  [Hedin,  1987]. 
However,  there  is  another  peak  just  before  dawn,  which  is  also  the  time  when  Hi  is  largest  (not 
shown).  The  pre-dawn  peak  is  indicative  of  the  fact  that  Ho  and  Hi  are  model  parameters,  taking 
on  large  magnitudes  at  a  time  when  the  electron  density  profile  broadens  as  a  result  of  low 
density.  This  low  density  is  due  to  the  cumulative  effect  of  recombination  in  the  bottomside  and 
diffusion  from  above,  which  results  in  a  broad-peaked  profile  shape. 

The  high  accuracy  of  the  linear-H  model  is  largely  due  to  the  fourth  parameter.  Hi.  This 
parameter  can  be  interpreted  as  the  average  altitude  gradient  in  the  pseudo-scale  height  above  the 
F2  peak  height.  The  parameter  Hi  distinguishes  the  linear-H  model  from  previous  formulations 
based  on  the  Chapman  profile  function  [Anderson  et  al.,  1987;  Burnside  et  al.,  1983].  Along  with 
Ho,  the  parameter  Hi  controls  the  shape  of  the  profile  above  the  F2  peak;  however.  Hi  =  0  below 
the  F2  peak.  The  units  of  Hi  are  km  per  km.  Figure  19  shows  that  the  distribution  of  logio  Hi  is 
skewed  negatively.  The  mode  of  the  logioHi  distribution  corresponds  to  Hi  of  about  6  km  per 
100  km  (Table  1).  The  parameter  Hi  peaks  at  about  8  km  per  100  km  around  08  AST,  and 
reaches  its  minimum  in  the  early  afternoon  (not  shown).  About  18%  of  the  profiles  have  Hi 
values  below  1  km  per  100  km.  The  mean  gradient  of  4  km  per  100  km  gives  the  necessary 
flexibility  for  the  linear-H  model  to  fit  the  Ne(h)  profiles  with  a  mean  accuracy  of  ±  5%  (Fig.  13, 
Table  1). 

The  fractional  uncertainties  in  Ho  and  Hi  are  distributed  as  shown  in  Figures  20  and  21, 
respectively.  The  uncertainties  are  obtained  from  the  diagonal  elements  of  the  model  covariance 
for  each  parameter,  using  equation  (13)  of  Picone  et  al.  (1997).  Two  thirds  of  the  a(Ho)/Ho 
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distribution  (Fig.  20)  fall  within  the  6%  level.  The  a(Ho)/Ho  di.stribution  peaks  at  between  2% 
and  3%  with  86%  of  the  fits  within  20%.  Uncertainties  larger  than  100%  have  a  cumulative 
frequency  of  less  than  0.13  (13%),  and  79%  of  the  Ho  uncertainties  are  10%  or  less.  The  mean 
becomes  5  ±  5%  if  only  profiles  with  uncertainties  within  100%  are  counted  (Table  1).  The 
percentage  uncertainty  in  Hj  (Fig.  21),  on  the  other  hand,  has  an  unusually  flat  distribution,  with  a 
peak  of  less  than  1%  at  a  rather  high  frequency  (Table  1).  In  contrast  to  Ho,  over  94%  of  the 
uncertainties  in  Hi  are  100%  or  less.  The  mean  decreases  to  19  ±  19%  if  the  latter  sub-population 
is  used.  Over  two  thirds  of  all  values  fall  within  the  23%  uncertainty  level.  There  is  a  greater 
range  of  uncertainties  in  a(Hi)/Hi  than  in  a(Ho)/Ho,  which  reveals  a  lower  sensitivity  of  the  data 
to  Hi  in  the  model.  However,  the  average  model  uncertainties  (Figures  13  and  14)  indicate  that 
the  overall  accuracy  of  the  model  is  better  than  that  shown  by  the  distributions  of  errors  in  Ho  and 
Hi  alone. 

3.4.  Distributions  in  the  Topside  Ionosphere  at  589  km 

There  is  considerable  interest  in  a  better  understanding  of  the  topside  region  above  the 
ionospheric  F2  layer.  Much  of  this  arises  from  the  need  to  improve  total  electron  content  (TEC) 
and  light  ion  models.  There  is  capability  to  detect  light  ion  fractions  with  the  Arecibo  incoherent 
scatter  radar.  The  reanalysis  of  the  1985-1995  World  days  applied  this  capability,  yielding  H"^ 
and  He'^  densities  (actually,  the  spectral  fits  provide  ion  fractions).  Since  the  standard  World  Day 
can  reach  as  high  as  680  km,  we  chose  the  range  gate  589  km  as  a  representative  altitude  at  which 
to  evaluate  the  measured  distributions  of  ion  fractions  as  well  as  electron  and  ion  temperatures. 
One  limitation  of  the  multi-frequency  technique  of  Sulzer  (1986)  is  that  the  bandwidth  is  reduced, 
compared  to  a  single-frequency,  in  order  to  sample  seven  simultaneous  spectra.  The  reduced 
bandwidth  can  limit  the  detection  of  large  H^  fractions  since  these  produce  frequency-broadened 
spectra.  The  broadened  spectra  can  bias  the  ion  temperature  determination  since  both  affect  the 
shape  of  the  spectrum  in  a  similar  way.  As  will  be  seen,  the  H^  fraction  distribution  at  589  km 
shows  relatively  low  fractions.  However,  some  of  the  measured  ratios  of  electron  and  ion 
temperature  (T/Tj)  are  less  than  unity,  suggesting  that  some  of  the  H*  fractions  in  the  CEDAR 
database  could  be  underestimating  the  actual  fractions.  Thus,  the  fractions  and  the  temperatures 
are  related  by  the  similarity  of  their  incoherent  scatter  radar  spectra. 
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We  present  the  first  long-term  distribution  of  Te  at  Arecibo,  shown  in  Figure  22.  The  F  region 
electron  and  ion  temperatures  are  important  in  determining  the  energy  budget,  the  electron  and 
ion  scale  heights,  and  in  understanding  the  overall  dynamic  behavior.  The  distribution  is  skewed 
towards  higher  temperatures,  exceeding  3000  K  with  a  cumulative  frequency  of  0.05.  However, 
99.1%  of  the  distribution  is  found  below  3500  K.  The  largest  and  smallest  temperatures  at  the 
0.01  level  are  3400  -  3500  K  and  500  -  600  K.  Thus,  Te  varies  by  as  much  as  a  factor  of  five  at 
the  1%  frequency  bin.  The  overall  low  frequency  levels  indicate  that  the  ionospheric  electrons  at 
this  altitude  occupy  a  wide  range  of  accessible  thermal  energies.  This  is  in  agreement  with  the 
significant  scatter  in  Te  found  by  Mahajan  (1967)  and  Kohnlein  (1981). 

The  electron  temperature  exhibits  a  bimodal  distribution  at  589  km,  with  two  broad  peaks 
separated  by  about  800  -  900  K,  beyond  the  range  of  the  standard  deviation.  This  bimodal 
distribution  is  due  to  two  factors.  The  first  is  the  diurnal  variation  in  which  nearly  all  of  the 
temperatures  lower  than  1500  K  are  nighttime  values.  The  second  is  that  Te  varies  with  changing 
solar  activity  [Grebowski  et  al.,  1993].  We  find  that  the  nighttime  Te  increases  with  solar 
activity  whereas  the  daytime  Te  decreases  with  increasing  solar  activity  (F10.7).  The  decrease 
with  increasing  F10.7  is  consistent  with  both  the  electron  density  at  589  km  and  the  peak  density 
increasing  with  higher  solar  activity.  In  summary,  the  topside  electrons  are  hotter  at  night  during 
high  solar  activity  but  are  cooler  during  daytime  whereas  the  low  FI 0.7  topside  is  colder  at  night 
but  hotter  during  daytime.  While  all  local  times  are  sampled  equally,  about  20%  of  the 
temperatures  during  low  solar  activity  were  measured  during  winter. 

The  ratios  of  electron  to  ion  temperature  at  589  km,  T/Ti,  are  shown  in  Figure  23.  There  is  a 
broad  maximum  between  1.0  and  1.3,  and  a  mean  of  1.3  ±  0.4;  the  median  is  1.22  (Table  1).  A 
ratio  of  one  would  be  expected  with  frequency  of  0.5  since  Tg  =  Ti  at  night,  and  nearly  half  of  all 
the  measurements  are  nighttime  (Fig.  4).  Carpenter  and  Bowhill,  [1971]  found  that  Tg/Tj  need 
not  be  one  throughout  the  entire  night,  particularly  during  winter.  Thus,  the  distribution  need  not 
be  evenly  distributed  about  Tg/Ti  =  1.  The  cumulative  frequency  or  population  for  Tg/Tj  >  1.1  is 
larger  than  0.6;  for  Tg/Tj  >1.3  the  frequency  is  about  0.4.  The  largest  Tg/Tj  ratios  occur  typically 
during  sunrise  when  the  plasma  is  least  dense  and  efficiently  heated  by  photoelectrons,  and  in 
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daytime  when  photoelectron  and  plasmaspheric  heating  rates  are  greatest.  At  589  km,  Tc/Tj  =  2  is 
seen  with  a  frequency  of  0.03  while  values  of  T/Tj  >  2  have  a  cumulative  frequency  of  close  to 
0.05.  However,  Tg/Ti  ratios  greater  than  two  can  be  regarded  as  outliers. 

Values  of  T/Ti  <  1  are  a  consequence  of  the  freedom  of  the  non-linear  least  squares  fit  of  the 
incoherent  scatter  spectra  to  assign  whatever  T/Ti  ratio  best  fits  the  measured  spectrum.  These 
low  ratios,  however,  are  unphysical  since  there  is  no  mechanism  allowing  the  more  massive  ions 
to  become  hotter  than  the  electrons.  Therefore,  T/Ti  <  1  values  are  an  artifact  of  the  fitting 
process  of  noisy  spectra,  though  as  many  as  13%  of  the  ratios  are  less  than  one.  Note  that  both 
day  and  night  values  of  T/Ti  are  shown  in  Fig.  23.  Fits  with  Tg/Tj  <  1  generally  are  found  at 
night  but  account  for  only  13%  of  the  distribution.  The  temperatures  and  the  and  He'^ 
fractions  are  related  in  the  radar  spectrum  [Evans  and  Louwenthal,  1964].  Low  levels  of  He"^  in 
the  spectrum  may  be  confused  for  high  ion  temperatures,  resulting  in  too  low  a  Tg/Tj  ratio. 
However,  a  cursory  examination  does  not  reveal  any  abnormality  in  the  light  ion  fractions 
correlating  with  Tg/T;  <  1.  In  addition,  the  spectral  reduction  assumes  all  ions  have  a  common 
temperature.  Recent  work  by  Sulzer  and  Gonzalez  (1996)  shows  distinct  and  H^  temperatures 
at  Arecibo,  and  these  improve  upon  the  spectral  fits.  Two-temperature  fits  were  not  done  with 
the  World  Day  dataset. 

The  distribution  of  the  measured  electron  density  at  589  km,  Ne(589  km),  is  shown  in  Figure  24. 

A  comparison  of  the  Ne(589  km)  and  nmax  distributions  (Fig.  6)  shows  that  the  former  has  a 
higher  (more  concentrated)  peak  than  the  nmax  distribution.  The  mean  in  Ne(589  km)  is  2.3  ± 

2.5  X  lO'  ‘  m'^,  with  frequency  of  about  0.01,  whereas  the  mean  of  the  nmax  distribution  is  9.8  ± 
6.7  X  lO' ' m'^,  with  similar  frequency  as  Ng  but  four  times  larger  bin  size  (Table  1).  A  smaller  bin 
size  decreases  all  the  frequencies,  resulting  in  a  less  frequent  peak  value.  The  Nc(589  km) 
distribution  is  more  peaked  than  nmax  in  spite  of  the  smaller  bin  size.  Another  way  to  examine 
the  spread  is  by  comparing  the  standard  deviation  as  a  fraction  of  the  mean  in  the  distribution: 

1.09  for  Ng  versus  0.68  for  nmax  (Table  1).  The  topside  density  is  more  spread  than  nmax.  The 
Ng  distribution  is  strongly  skewed  positive  though  there  are  only  156  values  larger  than  lO'^  m‘^ 
out  of  9801  sampled.  Thus,  there  are  very  few  outliers  but  a  greater  spread  of  electron  densities 
at  589  km  than  at  the  F2  peak  (Nmax).  Unlike  Tg  (Fig.  22),  however,  the  Ne  distribution  shows 
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no  bimodality.  This  suggests  that  the  bimodality  in  Te  does  not  require  a  corresponding 
bifurcation  in  the  electron  density.  The  skewness  in  Ne,  however,  is  large  enough  to  account  for 
the  bifurcation  in  Te(589  km). 

The  H'^  fraction  distribution  at  589  km,  shown  in  Figure  25,  is  weakly  bimodal  with  peaks 
separated  by  just  a  few  percent:  below  the  1%  fraction  and  just  above  2%.  The  distribution  is 
also  skewed  positive,  like  Ne(589  km)  (Fig.  24).  The  most  probable  density  can  be  estimated 
by  multiplying  the  modes  of  fraction  and  Ne  at  589  (Table  1)  or  about  0.75  x  10^  m'^  (750  cm' 
^).  The  individual  bin  frequencies  are  rather  low:  the  H^/Ne  =  6%  level  has  a  frequency  of  only 
0.02.  The  mean,  <H^/Ne>  =  0. 13  +  0.21,  has  a  frequency  of  only  0.02,  and  the  median  fraction  is 
just  0.04.  Proton  fractions  of  15%  and  higher  have  a  cumulative  frequency  of  0.23  whereas 
fractions  greater  than  50%  (i.e.,  0.5)  have  a  cumulative  frequency  of  0.07.  Overall,  fractions 
between  the  0.10  and  0.25  fractions  are  about  equally  likely.  Higher  proton  fractions  are  much 
less  frequent  in  the  distribution.  The  skew  in  the  distribution  is  due  to  the  variability  of  the 
density  in  the  topside,  which  varies  with  local  time  [Tepley  and  Kerr,  1987;  Ho  and  Moorcroft, 
1971],  magnetic  activity  [Tepley  and  Kerr,  1987],  season  [Tepley  and  Kerr,  1987;  Prasad,  1968], 
and,  solar  cycle  [Bailey  et  al.,  1982;  Evans  and  Holt,  1978;  Prasad,  1968]. 

Modeling  [Bailey  et  al.,  1982]  shows  that  the  H^  fraction  is  most  sensitive  to  changes  in  the 
neutral  hydrogen  density,  which  itself  varies  with  solar  cycle  by  as  much  as  a  factor  of  two 
[Tinsley,  1978].  At  a  fixed  altitude,  H^/Ne  increases  during  solar  minimum  since  neutral 
hydrogen  is  more  abundant  than  during  solar  maximum  conditions  at  all  latitudes  [Kutiev  et  al., 
1980].  Since  H"^  is  produced  primarily  by  charge  exchange  between  and  H,  the  large 
abundance  of  hydrogen  in  solar  minimum  results  in  a  greater  proton  fraction.  However,  diurnal 
changes  in  neutral  hydrogen  can  be  just  as  large.  Thus,  the  flat  portion  of  the  H'^/Ne  distribution 
contains  both  solar  cycle  and  diurnal  variations.  In  addition,  the  cumulative  frequency  of  the  50% 
fraction  or  higher  is  of  only  0.07.  This  low  frequency  suggests  that  the  O"^— H^  transition  height  is 
rarely  found  at  589  km  or  that  the  multi-frequency  technique  is  limited  in  its  ability  to  sample 
high  proton  fractions.  The  transition  height  can  be  as  low  as  475  km  at  night  at  Arecibo 
[Moorcroft,  1969;  Prasad,  1968;  Carlson  and  Gordon,  1966].  However,  such  low  transition 
heights  are  generally  found  in  the  early  morning  near  0400  AST  during  solar  minimum  winter. 
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Also,  OGO  6  [Kutiev  et  al.,  1980]  and  AE  E  satellite  [Gonzalez  et  al.,  1992]  measurements,  as 
well  as  constrained  incoherent  scatter  radar  measurements  at  Arecibo  [Erickson  and  Swartz, 
1994],  find  the  transition  height  generally  above  600  km. 

The  He'^/Ne  distribution,  shown  in  Figure  26,  has  a  large  mode  between  0%  and  0.5%,  unlike 
H^/Ne.  Low  He"^  fractions  are  expected  from  other  observations  and  modeling  [Bauer,  1966; 
Prasad,  1968;  Moorcroft,  1969;  Bailey  and  Sellek,  1990;  Hoegy  et  al.,  1991;  Gonzalez  et  al., 
1992],  and  from  the  fact  that  neutral  helium  is  less  abundant  than  neutral  hydrogen  in  the 
thermosphere  and  exosphere.  Even  the  recently  discovered  He"^  bulge  at  Arecibo  does  not  appear 
to  exceed  the  25%  fraction  [Gonzalez  and  Sulzer,  1996].  The  mean  of  the  entire  distribution, 
<He^/Ne>  =  0.03,  occurs  with  frequency  of  less  than  0.05;  the  standard  deviation  of  0.09  is  three 
times  the  size  of  the  mean.  The  He"^  distribution  is  more  narrowly  peaked  than  that  of  H^,  so  the 
two  are  quite  different.  Besides  the  mode,  the  next  bin  is  about  10  times  less  frequent  than  the 
mode.  This  is  due  to  the  large  number  of  “zero”  fractions  in  the  distributions  (the  median  is 
0.0000).  A  zero  fraction  does  not  mean  zero  He"^;  instead,  the  actual  fraction  is  too  low  or  too 
high  to  be  resolved  in  the  radar  spectral  fits.  The  ion  fraction  distributions  at  589  km  show  that; 
(i)  He‘*^/Ne  is  generally  smaller  than  H^/Ne;  (ii)  both  are  skewed  positive;  and,  (iii)  He'^/Ne  has 
lower  frequencies  beyond  10%. 

The  difference  in  the  two  distributions  can  be  explained  by  the  difference  in  the  production  and 
loss  processes.  Production  of  He^  is  largely  due  to  photoionization  (by  solar  EUV  photons  below 
50.4  nm)  at  high  altitudes  whereas  the  loss  rate  peaks  at  low  altitudes.  Above  the  F2  layer,  the 
loss  rate  is  small  since  there  is  little  N2  for  He^  to  recombine.  The  largest  He^  fraction  and 
density  occurs  at  night  due  to  inflow  of  He^  from  the  plasmasphere.  Modeling  by  Bailey  and 
Sellek  (1990)  shows  that  HeV  >  1  at  500  km  prior  to  sunrise  in  the  L=1.8  flux  tube  during 
solar  minimum  equinox.  The  distributions  here  obtained  provide  some  support  for  the  possibility 
of  equal  fractions  and  concentrations:  the  10%  fraction  has  roughly  equal  probability  for  both 
ions.  Of  course,  that  does  not  establish  coincidence,  and,  in  an  aggregate  view,  large  He^ 
fractions  have  lower  frequencies  than  the  corresponding  fractions.  The  experimental 
uncertainties  in  the  two  ions  are  not  equal  for  the  same  fractional  values  obtained  from  fitting  the 
radar  spectra  (the  He^/Ne  uncertainties  are  generally  higher).  Another  difference,  discussed 
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previously,  is  the  small  negative  bias  present  in  the  He^  distribution,  which  is  an  artifact  of  the 
fitting  process  that  occurs  with  very  noisy  spectra  when  a  common  ion  temperature  is  assumed. 
The  bias  is  removed  by  allowing  individual  ions  to  have  separate  temperatures,  as  done  by  Sulzer 
and  Gonzalez  [1996].  There  is  not  enough  frequency  bandwidth  in  the  World  day  spectra  to  do 
multiple  ion  temperature  fits. 

The  Arecibo  World  Days  measured  the  radial  (line-of-sight)  ion  velocity  component  as  a  function 
of  altitude  and  beam  azimuth.  The  radial  speeds  are  used  in  obtaining  the  vector  ion  velocity  and 
the  ExB-drifts  of  the  plasma  [Behnke  and  Haiper,  1973;  Sulzer,  1986;  Burnside  et  al.,  1987; 
Fejer,  1993].  A  single  altitudinal  profile  of  radial  speeds  is  obtained  at  each  azimuth  position  for 
a  total  of  12  profiles  every  360°  swing  of  the  radar  beam.  The  radial  speeds  at  589  km  have  a 
broad  but  symmetric  distribution,  shown  in  Figure  27,  centered  at  zero  velocity,  with  mean  of  -4 
m  s'*  (Table  1).  The  mean  becomes  -0.09  ±  31  m  s'*  if  only  velocities  within  ±  100  m  s'*  are 
considered,  which  comprise  95%  of  the  population.  The  median  in  either  case  is  zero  velocity. 
There  is  no  net  plasma  flow  across  589  km  on  average,  and  positive  and  negative  velocities  are 
about  equally  probable.  The  large  variability  is  evidenced  by  the  broad  wings  in  the  distribution. 
There  are  certainly  periods  of  net  plasma  transport  in  or  out  (down  or  up)  but  these  are  balanced 
in  the  long-term. 

That  the  radial  speed  is  a  measure  of  net  plasma  motion  across  the  589  km  boundary  is  both 
intuitive  and  formal.  The  line-of-sight  velocity  at  height  z,  azimuth  <j)  (measured  counter 
clockwise  from  magnetic  east),  and  zenith  angle  0,  Vr,  is  given  by 

Vr(z,  0,  (}))  =  Sin0  [  Vs  Sin<()  +  Ve  COSij)  ]  +  Vz  COS0  (7) 

where  Vj  is  the  horizontal  ion  velocity  along  the  magnetic  southward  direction,  Ve  is  the  eastward 
ExB  drift  velocity  component,  and  v^  is  the  vertical  ion  drift  component  due  to  vertical  winds  and 
given  by 

Vz  =  Uz  sin  I  (8) 
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where  is  the  vertical  wind  component,  and  I  is  the  magnetic  inclination  angle.  The  southward 
ion  drift,  Vs,  is  given  by 

Vs  =  vl  cos  I  -  VxN  sin  I  (9) 

where  vl  is  the  parallel  ion  velocity  (positive  up  along  the  magnetic  field  line),  and  v^n  is  the 
ExB-drift  in  the  perpendicular-north  direction.  The  parallel  ion  velocity  contains  the  influence  of 
ambipolar  diffusion  [Schunk  and  Walker,  1970],  which,  along  with  meridional  neutral  winds, 
drive  plasma  motions  along  the  field  line.  Plasma  motions  along  the  radar  line-of-sight  are  the 
sum  of  these  velocities  which,  with  the  possible  exception  of  vertical  neutral  winds,  are  present 
with  varying  magnitudes  [Burnside  et  al.,  1981].  A  downward  radial  velocity  arises  from  either  a 
downward  parallel  drift,  or  a  downward  ExB-drift  component,  or  a  linear  combination.  Zonal 
ExB-drifts  may  cause  radial  velocities  since  there  is  a  component  along  the  line-of-sight 
whenever  the  beam  points  away  from  the  magnetic  meridian.  The  zonal  drift,  however,  does  not 
move  plasma  vertically.  Zonal  drifts  will  bias  the  radial  distribution  only  if  there  is  a  non-zero 
zonal  gradient  in  the  zonal  ExB  drift.  A  bias  can  be  also  caused  by  a  substantial  change  in  the 
electric  field  during  the  beam  swing  (16-min). 

The  fact  that  the  distribution  is  symmetrical  about  zero  with  near  zero  median  implies  that  the 
zonal  gradient  is  near  zero  or  constant  on  average.  According  to  Equation  (8),  a  non-zero  but 
constant  zonal  gradient  induces  a  radial  velocity  when  the  beam  points  east  but  an  equal  and 
opposite  component  occurs  with  the  beam  in  the  western  azimuth  sector.  The  resulting 
components  contribute  equal  positive  and  negative  velocities  but  do  not  bias  the  azimuthal 
distribution.  Only  actual  vertical  velocities  can  generate  a  bias  or  moment  in  the  distribution. 
Therefore,  the  zero  bias  cannot  be  caused  by  zonal  gradients  in  the  zonal  drift,  and  the 
distribution  of  Vr  is  evidence  of  no  net  vertical  plasma  transport  across  589  km  over  long  time 
averages. 

3.5.  Ap  Index 


20 


The  level  of  geomagnetic  activity  in  the  dataset  can  be  seen  in  the  distribution  of  the  Ap  index, 
shown  in  Figure  28.  The  Ap  index  is  a  global  average  of  the  eight  daily  values  of  ap,  which  is  the 
K  index  converted  to  a  linear  scale.  The  average  is  made  over  all  the  stations  reporting  K  values, 
located  mostly  in  the  northern  hemisphere.  The  Ap  index  ranges  between  0  and  400;  Ap  <  15 
(Kp  <  3)  denotes  generally  quiet  conditions,  whereas  Ap  >  48  (Kp  >  5)  denotes  disturbed 
conditions.  Most  of  the  data  has  Ap  <  20;  the  mode  is  6,  and  the  median  is  1 1.  The  mean  Ap  is 
14  ±  12,  with  frequency  of  0.06.  Values  of  Ap  >  27  (Kp  =  4)  have  a  cumulative  frequency  of 
less  than  0.04.  Therefore,  less  than  4%  of  all  the  World  days  occurred  during  Ap  values 
corresponding  to  geomagnetically  disturbed  conditions.  Disturbed  conditions  may  account  for 
4%  frequency  level  effects  in  the  distributions. 


4.  Summary  and  Conclusions 

We  present  an  “ionospheric  climatology”  comprised  of  histograms  of  10  years  of  World  Day 
measurements  obtained  with  the  Arecibo  incoherent  scatter  radar  as  stored  in  the  NCAR-CEDAR 
database.  In  addition,  we  compare  the  linear-H  and  splines  model  fits  to  the  measured  electron 
density  profiles,  the  latter  model  being  the  one  used  in  obtaining  the  peak  height  and  densities  in 
the  CEDAR  database. 

1.  About  55%  of  the  observations  correspond  to  solar  10.7-cm  flux  of  120  and  higher.  All  local 
times  were  sampled  equally  (to  within  20%)  in  the  distribution  of  about  140,000  profiles 
obtained  by  pointing  the  beam  at  an  elevation  angle  of  75°. 

2.  Model  fits  to  the  measured  electron  density  profiles  in  the  database  were  made  using  the  four- 
parameter  linear-H  model  consisting  of  F2  peak  density  (nmax),  height  (hmax),  scale  height 
(Ho)  and  linear  scale  height  gradient  (Hi).  As  many  as  1 1,440  profiles  were  fitted,  90%  with 
better  than  10%  accuracy  on  average.  Thus,  the  linear-H  model  is  accurate  and  precise  to 
within  10%.  The  model  accuracy  in  fitting  measured  electron  density  profiles  matches  the 
experimental  uncertainty  of  10%  at  the  “one-sigma”  level  or  higher.  The  four-parameters  of 
the  linear-H  model  retain  physical  meaning,  although  the  fitting  procedure  treats  them  as 
purely  numerical  parameters. 
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3.  Both  the  model  and  measured  (cubic-splines)  peak  density  distributions  are  very  similar. 
Overall,  the  linear-H  model  peak  densities  below  about  0.3  x  lO'^  m'^  are  lower  than  those  in 
the  CEDAR  database.  Above  this  values,  the  linear-H  peak  densities  are  up  to  20%  to  30% 
higher,  as  seen  by  comparing  the  mean  and  medians  (Table  1)  of  the  distributions  of  the  peak 
densities  from  linear-H  and  from  the  CEDAR  database  (Fig.  6).  A  more  precise  measure  of 
this  bias  is  seen  in  the  distribution  of  residuals  between  the  linear-H  and  database  peak 
densities,  which  shows  that  negative  residuals  are  more  frequent  than  positive  residuals  (Fig. 
11). 

4.  The  linear-H  distribution  of  peak  heights  is  highly  Gaussian.  The  distribution  of  peak 
heights  obtained  with  splines,  on  the  other  hand,  “fluctuates”  about  the  corresponding  linear- 
H  distribution  (Fig.  7),  as  expected  from  the  fact  that  the  splines  do  not  account  for 
measurement  noise.  The  mean  in  the  CEDAR  peak  height  (zmax)  is  320  ±  49  km  whereas 
the  linear-H  peak  height  (hmax)  is  316  ±  51  km;  the  corresponding  medians  are  317  and  316 
km,  respectively  (Table  1).  Therefore,  the  linear-H  peak  height  is  a  few  km  lower  than  the 
splines  zmax.  In  fact,  the  model  underestimates  zmax  by  about  1  km  on  average  as  seen  from 
the  mean  and  median  of  the  distribution  of  residuals  between  the  linear-H  and  CEDAR 
database  values  (Fig.  12).  Nearly  68%  of  that  distribution  has  negative  residuals,  confirming 
that  the  linear-H  hmax  is  smaller  than  the  CEDAR  database  value,  obtained  using  a  splines 
fit.  In  general,  the  linear-H  model  provides  slightly  lower  peak  height  and  densities  than  a 
splines  fit. 

5.  The  (reduced)  of  the  linear-H  fits  to  the  measured  electron  density  profiles  are  distributed 
about  one,  which  is  the  value  expected  of  optimal  fits.  The  distribution  of  logioX^  is  a  broad 
Gaussian  with  standard  deviation  of  ±  6  about  a  mean  of  1.15,  with  84%  of  the  fits  having  x^ 
values  between  0.1  and  10.  A  broad  distribution  can  result  from  a  p.seudo-random  or  partial 
bias  in  the  uncertainties  assigned  to  the  data. 

6.  The  model  precision  in  both  peak  density  and  peak  height  is  within  10%  in  about  85%  of  the 
fits.  The  linear  uncertainties  in  hmax  are  more  narrowly  distributed  than  nmax,  as  seen  by 
comparing  Figures  9  and  10.  The  fact  that  these  errors  have  normal  distributions  preclude  any 
significant  model  bias  in  the  calculated  variance. 
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7.  The  linear-H  peak  density  (nmax)  is  within  10%  of  the  CEDAR  database  values  (Nmax)  in 
over  96%  of  the  fits,  and  within  5%  in  over  80%  of  the  fits.  The  model  underestimates  the 
Nmax  values  obtained  with  splines  by  less  than  1%  on  average,  as  seen  by  the  distribution  of 
the  residuals  between  the  model  and  measured  peak  densities  (Fig.  1 1).  Negative  residuals 
are  cumulatively  more  frequent,  as  seen  from  median  and  mean  values  of  nearly  -1%,  and 
from  the  negative  skew  in  the  distribution. 

8.  The  linear-H  hmax  is  within  10%  of  the  CEDAR  values  in  about  98%  of  the  fits,  and  within 
5%  in  over  90%  of  the  fits.  The  distribution  of  (normalized)  peak  height  residuals  (Fig.  12)  is 
narrower  than  the  corresponding  peak  density  residuals  (Fig.  11).  The  median  of  the 
distribution  of  residual  is  -1.2%,  a  slight  underestimation  by  the  linear-H  model.  The  linear-H 
peak  height  and  densities  have  similar  accuracy,  with  virtually  all  residuals  within  ±10%. 

9.  The  overall  accuracy  of  the  fits  to  the  measured  electron  density  profiles  is  evaluated  from  the 
distribution  of  the  average  absolute  residual  for  each  profile,  as  a  fraction  of  the  measured 
density  (Fig.  13).  By  this  measure,  the  linear-H  profiles  are  well  within  10%  of  the  data  in 
94%  of  the  fits,  and  within  5%  in  66%  of  the  fits.  The  linear-H  model,  therefore,  fits  the 
measured  electron  density  profile  with  high  accuracy  (10%  or  better).  The  model  accuracy  is 
comparable  to  the  average  precision  of  the  measured  density  profile  (Fig.  14),  which  is  within 
10%  in  89%  of  the  profiles,  and  within  5%  in  56%  of  the  profiles.  The  median  experimental 
precision  is  larger  (4.6%)  than  linear-H  (3.6%)  by  only  1%.  Therefore,  the  model  and  the  data 
have  nearly  identical  accuracy  and  precision.  Furthermore,  the  distribution  of  the  difference 
between  the  linear-H  residuals  and  the  experimental  uncertainties,  also  the  difference  of  the 
previous  two  errors,  shows  that  94%  of  all  the  fits  are  within  10%  of  the  experimental 
precision,  and  71%  of  the  profiles  are  within  5%.  Therefore,  the  average  linear-H  residual  is 
within  10%  of  the  measured  profile  error. 

10.  The  linear-H  model  parameter  Ho  takes  on  values  consistent  with  those  expected  from  the 
atomic  oxygen  scale  height  near  the  F2  peak,  namely  a  mean,  median  and  mode  of  about  50 
km  with  a  standard  deviation  of  16  km.  The  Hq  distribution  can  be  regarded  as  a  proxy  for  the 
neutral  temperature  near  the  peak,  since  the  atomic  oxygen  scale  height  is  directly 
proportional  to  the  temperature.  The  most  frequent  uncertainty  in  Ho,  as  well  as  the  median 
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uncertainty,  is  4%;  also,  79%  of  the  uncertainties  are  within  10%.  Furthermore,  low  values 
of  Ho  of  20  km  occur  mostly  for  F10.7  values  of  100  or  less.  The  distribution  of  the  preeision 
in  Ho  is  similar  to  the  average  electron  density  profile  error  (Figs.  20  and  14). 

11.  The  fourth  parameter  in  the  linear-H  model.  Hi,  is  the  average  linear  altitude  gradient,  which 
modifies  the  pseudo-scale  height  Ho  and  the  electron  density  profile  shape  above  the  F2  layer 
height.  About  82%  of  all  the  model  profiles  have  Hi  values  greater  than  0.01  scale  heights 
per  km  (Fig.  19).  The  average  gradient  in  the  model  scale  height  is  6  km  /  100  km  (Table  1). 
Therefore,  the  linear-H  model  implies  non-zero  scale  height,  and  temperature,  variation  with 
altitude  above  the  F2  peak.  The  mode  in  the  distribution  of  the  precision  in  H|  is  less  than  1% 
though  the  distribution  is  strongly  skewed  positive,  with  66%  of  the  errors  within  23%.  Over 
94%  of  the  errors  are  of  100%  or  less.  The  model  assigns  large  errors  to  Hi  even  when  the 
overall  fit  uncertainties  are  small.  Moreover,  the  skew  in  the  Hi  precision  contrasts  with  the 
more  symmetrical  distribution  of  the  precision  in  Hq  (Fig.  20). 

12.  The  electron  temperature  at  589  km  has  a  strong  bimodal  distribution,  which  is  consistent 
with  solar  cycle  and  local  time  dependence.  The  histogram  peaks  are  separated  by  about  900 
K  and  peak-to-valley  frequency  ratios  of  almost  a  factor  of  three.  The  database  shows  99%  of 
the  Te  distribution  falling  below  3500  K.  The  Te/Tj  ratio  at  589  km  (Fig.  23),  on  the  other 
hand,  has  a  unimodal  distribution  with  a  peak  at  Te/Ti  =1.03.  A  single  Te/Ti  peak  suggests 
that  the  ion  temperature  is  also  bimodal,  following  the  electron  temperature  closely.  Values 
of  Te/Ti  <  1  appear  in  13%  of  the  distribution  but  these  are  unphysical,  and  due  to  the  freedom 
afforded  to  the  radar  spectral  fits.  The  topside  electrons  are  hotter  at  night  during  high  solar 
activity  but  are  cooler  during  daytime,  whereas  the  low  F10.7  topside  is  colder  at  night  but 
hotter  during  daytime. 

13.  The  electron  density  at  589  km  is  unimodal,  unlike  the  electron  temperature.  The  distribution 
is  flat  between  0.4  x  lO'^  and  1  x  lO'^  m  '\  The  density  at  589  km  has  a  stronger  (more 
frequent)  mode  than  the  peak  density  (Fig.  6),  which  suggests  a  different  variability  occurs  in 
the  topside  electron  density.  On  the  other  hand,  the  spread  of  densities  at  589  km  is  larger 
than  nmax.  As  a  fraction  of  the  mean,  the  standard  deviation  at  589  km  is  about  60%  larger 
than  the  standard  deviation  in  nmax.  Therefore,  there  is  a  wider  spread  of  topside  densities  in 
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Ne  than  in  nmax  though  the  topside  distribution  is  more  sharply  peaked.  The  bimodality  in 
the  electron  temperature  (Fig.  22)  is  not  seen  in  the  electron  density. 

14.  The  light  ion  (H^  and  He^)  fractions  at  589  km  are  skewed  positive.  The  proton  fraction 
distribution  is  weakly  bimodal,  with  a  peak  of  less  than  1%  and  another  at  2.5%.  The  most 
probable  density  is  7.5  x  lO’ m'^  (750  cm'^).  Proton  fractions  of  0.25  (25%)  and  higher 
occur  in  16%  of  the  data  while  fractions  between  10%  and  25%  have  about  equal  probability. 
No  significant  correlation  is  found  between  T/Tj  less  than  one  and  particular  or  He^ 
fractions.  The  distribution  of  He^  fractions  has  a  strong  peak  at  about  zero  with  frequency  of 
over  0.6,  which  is  narrower  than  the  proton  fraction.  The  wider  He"^  distribution  is  also  seen 
by  the  fact  that  the  standard  deviation  is  three  times  larger  than  the  mean,  whereas  the 
standard  deviation  is  only  1.61  larger  than  its  mean.  The  median  of  the  He'^  fraction 
distribution  is  0.00000,  since  half  of  the  spectral  fits  did  not  find  He"^.  The  mean  He^/Ne  is 
0.03  versus  0.13  for  H^/Ne,  that  is,  the  radar  sees  four  times  more  H^than  He^.  Thus,  the  two 
ion  fractions  have  different  distributions,  consistent  with  having  different  production  and  loss 
processes. 

15.  The  distribution  of  radial  (line-of-sight)  velocities  at  589  km  is  symmetric  about  zero,  with 
95%  of  the  velocities  within  ±  100  m  s’'.  The  median  velocity  is  zero,  and  the  mean  is  -0.09 
m  s'*.  Therefore,  there  is  no  net  plasma  flow  across  589  km  on  average.  Zero  velocity  is  the 
most  frequent  value. 

16.,  About  96%  of  all  the  data  corresponds  to  Ap  values  of  27  of  less,  corresponding  to  quiet 
conditions  (Kp  <  4).  Therefore,  only  distribution  values  with  cumulative  frequencies  of  less 
than  4%  can  be  directly  attributed  to  high  geomagnetic  activity. 
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Figures 


1.  Distribution  of  daily  I0.7-cm  solar  flux  (F10.7)  values  corresponding  to  the  dataset  of  ten  years  of  Arecibo 

World  days  between  1985  and  1995,  About  55%  of  the  distribution  has  FI 0.7  greater  than  120. 

2.  Distribution  of  years  covered  in  the  dataset  beginning  late  in  1985  and  extending  to  1995. 

3.  Distribution  of  months  covered  in  the  Arecibo  World  day  database  between  1985  and  1995.  January  is  month 
number  one. 

4.  Distribution  of  local  times  corresponding  to  the  fits  to  ten  years  of  Arecibo  World  days  between  1985  and 
1995.  Midnight  is  0.  Nearly  all  local  times  were  sampled  equally. 

5.  Virtually  all  ten  years  of  World  day  observations  were  made  with  the  antenna  beam  pointed  at  an  elevation 
angle  of  75°  (zenith  angle  of  15°)  while  the  azimuth  angle  was  changed. 

6.  Two  distributions  of  Fp  peak  densities  are  shown:  the  linear-H  fits  (solid  line,  “nmax”)  to  the  measured 
electron  density  profiles  contained  in  the  database,  and  splines  fits  (broken  line,  “Nmax”).  The  horizontal 
scale  is  in  electrons  per  10^  cm'^  (10*^  m'^).  The  CEDAR  database  Nmax  are  obtained  with  cubic  splines. 

7.  Two  distributions  of  F2  peak  heights  are  shown:  the  linear-H  fits  (solid  line,  “hmax”)  to  the  measured  electron 
density  profiles  contained  in  the  CEDAR  database,  and  the  splines  fits  (dotted  line,  '‘zmax”).  The  horizontal 
scale  is  the  peak  height  in  km.  Splines  do  not  account  for  random  experimental  noise. 

8.  Distribution  of  the  logio  (reduced  x^)  values  obtained  in  fitting  9801  electron  density  profiles.  The  distribution 
is  centered  at  0  (x^  =  1 ). 

9.  Distribution  of  the  percentage  model  uncertainty  in  the  linear-H  F2  peak  density  (nmax). 

10.  Distribution  of  the  percentage  uncertainty  in  the  model  F2  layer  height  (hmax). 

1 1 .  Distribution  of  the  difference  between  the  linear-H  and  CEDAR  database  peak  densities,  normalized  by 

Nmax  (CEDAR  database,  splines),  expressed  in  percent. 

12.  Distribution  of  the  difference  between  the  linear-H  and  CEDAR  database  peak  heights,  normalized  by  zmax 
(CEDAR  database,  splines),  expressed  in  percent. 

13.  Distribution  of  B,  the  absolute  magnitude  of  the  average  residual  between  the  linear-H  and  measured  electron 
densities  for  each  of  the  9801  fitted  profiles,  normalized  by  the  measured  electron  densities  (Equation  4), 
expressed  as  percent.  This  parameter  is  a  measure  of  the  accuracy  and  bias  of  the  linear-H  fit  to  the  measured 
profile. 

14.  Distribution  of  D,  the  average  of  the  measurement  uncertainties  in  each  electron  density  profile,  normalized  by 
the  data,  and  expressed  as  percent  (Equation  5).  This  is  a  measure  of  the  range  of  data  errors,  closely 
resembling  the  distribution  of  Fig.  13. 

15.  The  distribution  of  F  -  D,  the  difference  between  the  average  model  and  measured  residuals  for  each  profile 
(Fig.  13)  and  the  average  measurement  error  in  a  profile  (Fig.  14).  Negative  values  indicate  that  experimental 
errors  dominate. 

16.  Distribution  of  F,  the  median  of,  the  absolute  magnitude  of  the  residual  between  the  linear-H  and  measured 
electron  density  profiles,  normalized  by  the  measurement,  and  expressed  as  percent  (Equation  6).  This  is  a 
measure  of  the  median  fit  residual,  which  is  less  sensitive  to  outliers. 
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17.  Distribution  of  C,  the  minima  in  the  ratio  of  the  experimental  uncertainty  to  the  !inear-H  electron  density  for 

each  profile  (Equation  7).  This  is  a  measure  of  the  smallest  experimental  uncertainty  as  a  percentage  of  the 

linear-H  density. 

18.  Distribution  of  the  linear-H  parameter  Ho  (km),  a  pseudo-scale  height  term,  which  controls  the  shape  of  the 
model  electron  density  profile  (Equation  1).  This  parameter  takes  on  values  resembling  the  atomic  oxygen 
scale  height. 

19.  Distribution  of  the  linear-H  model  parameter,  logio  Hi,  which  modifies  Hq  and  thus  the  shape  of  the  linear-H 
electron  density  profile  above  hmax  (Equation  1).  The  parameter  Hi  is  the  linear  gradient  in  the  pseudo-scale 
height,  which  assumes  a  particular  value  for  each  model  profile. 

20.  Distribution  of  the  percentage  error  in  Ho  (the  linear-H  uncertainty  divided  by  Hq). 

21.  Distribution  of  the  percentage  error  in  Hi  (the  linear-H  uncertainty  divided  by  Hi). 

22.  Distribution  of  the  electron  temperature  measured  at  589  km  as  recorded  in  the  dataset  of  ten  years  of  Arecibo 
World  days.  The  horizontal  scale  is  degrees  K.  The  Tg  distribution  is  bimodal,  a  feature  consistent  with  solar 
cycle  and  diurnal  control. 

23.  Distribution  of  the  ratio  of  electron  to  ion  temperatures,  T/Tj,  measured  at  589  km.  Values  of  T/Tj  less  than 
one  comprise  are  an  artifact  of  the  radar  spectral  fits.  Only  13%  of  the  ratios  are  smaller  than  one. 

24.  Distribution  of  electron  densities  measured  at  589  km.  The  horizontal  scale  is  in  units  of  electrons  per  10*^ 
m‘^  (10^  cm’^). 

25.  Distribution  of  proton  fractions  (H'^/Ne)  measured  at  589  km. 

26.  Distribution  of  Helium  ion  fractions  (He'^/Ne)  measured  at  589  km. 

27.  Distribution  of  radial  (line-of-sight)  ion  velocities  along  the  antenna  beam  as  measured  at  589  km  (see 
equation  7).  The  distribution  is  centered  at  zero  velocity,  indicating  no  net  up  or  down  ion  motion.  Positive 
velocities  are  up  along  the  beam 

28.  Distribution  of  Ap,  the  global  average  of  the  3-hourly  ap  values.  Ap  values  of  27  correspond  to  Kp  =  4;  Ap  = 
48  corresponds  to  Kp  =  5,  which  denotes  magnetically  disturbed  conditions.  About  4%  of  the  World  day 
profiles  correspond  to  Ap  values  indicating  disturbed  conditions. 

29.  Composite  panel  of  all  the  figures  previously  shown. 


31 


year 

Figure —  2 


33 


DMA\pcuape\\f:\rdp8vr\arecibo_code\PLOT-AREClBO^OKjr.PRO\\(ONE-HISTOGRAM)\\ie-Nov-1098  15:63 


Fraction  of 
population  in  bin 


Elevation 
Figure —  5 


36 


DMA\pcuap6\\f:\rdpsvr\areclbo.code\PLOT_ARECIBO.OK-.T.PRO\\{ONELHISTOGRAM)\\ie-Nov-1998  15:53 


Fraction  of 
population  in  bin 


0.040 


0.030 


0.020 


0.010 


0.000 

li 


■8 


DMA\pcuap6\\f:\rdp3vr\arecIbo.^otie\PLOT-ARECIBO_OK.T.PRO\\(ONE-J«STOGRAM)\\l6-Nov-l998  15:53 


Fraction  of 
population  in  bin 


Figure —  9 


40 


DMA\pcuap6\\f;\rdp8vr\arcclbo_code\PL0T^RECIB0^0ICT.PR0\\(0NE_HIST0GRAM|\\l6-Nov-1998  15:53 


0  1  2  3  4  5  6 

sigma(hmax)/hmax,  % 

Figure —  10 


41 


DMA\peuap0\\f:\rdpsvr\areclbo„code\PLOT^RECIBO_OKJr.PRO\\(ONE^ISTOGRAM)\\16-Nov-1998  16:63 


Fraction  of 
population  in  bin 


42 


DMA\pcuap6\\f:\rdpsvr\arecibo^code\PLOT_ARECIBO>OK_T.PRO\\(ONE-HISTOGRAM)\\ie-Nov-1990  15:53 


Fraction  of 
population  in  bin 


-10  -5  0  5  10  15 

hmax  (Fit-CEDAR)/(CEDAR),  % 

Figure —  12 


43 


DMA\pcuap6\\f:\rdpsvr\areclbo_code\PL0T_ARECIB0>0ICT.PR0\\(0NE-HIST0GRAM)\\16-Nov- 1998  1 5:63 


Fraction  of 
population  in  bin 


I  I  I 

25 


44 


DMA\pcuap6\\f:\rdpsvr\arecIbo_code\PL0T_AREClB0_0K_T.PR0\\(0NE_HIST0GRAM)\\l6-Nov-1998  15:53 


Fraction  of 
population  in  bin 


46 


DMA\pcuap6\\f:\rdpavr\arecibo«code\PL0T^REClB0.0K_T.PR0\\(0NELHIST0GRAM)\\16-Nov7l998  15:53 


Fraction  of 
population  in  bin 


median  fit  error,  % 
Figure —  16 


47 


DMA\peuap6\\f:\rdpsvr\arecibo_eode\PL0T-ARECIB0.0{CT.PR0\\(0NE-HIST0GRAM)\\l6- Nov- 1990  16:64 


Fraction  of 
population  in  bin 


49 


DMA\pcuap6\\f;\r<ipavr\arecibo_code\PL0T^RECIB0^0KJr.PR0\\{0NE_HIST0CRAM)\\l8-Nov-1998  15:64 


Fraction  of 
population  in  bin 


-2.0  -1.5  -1.0  -0.5  0.0 

Fit  log(H,) 

Figure —  19 


j _ I _ I 

0.5 


50 


DMA\pcuap6Wf;\rdpsvr\arecibo_code\PLOT-ARECtBO_OK_T.PROW(ONE_HlSTOGRAM)'y^l6-Nov-199B  15:54 


Fraction  of 
population  in  bin 


51 


DMA\peuftp6\\f:\rdp8vr\ttPeclbO-eode\PL0T^RECIB0-0K_T.PR0\\(0NE_HIST0GRAM)\\16-Nov-1998  16:54 


Ne(  589  km  ),  10*  cm  * 
Figure —  24 


55 


DMA\pcuap0\\f:\rdp8vr\areclbo-.code\PLOT_ARECIBO-OKJr.PRO\\(ONE_HISTOGRAM)\\16-Nov-l998  15:54 


-0.1  0.0  0.1  0.2  0.3 

HeVNe,  589  km 

Figure —  26 


57 


DMA\pcuap6\\f:\rdp8vr\areclbo-,code\PLOT-AREC!BO-OIC.T.PRO\\(ONE_HISTOGRAM)\\ie-Nov-ie90  15:64 


Fraction  of 
population  in  bin 


59 


DMA\peuap6W:\rdp*vr\arecibo_eode\PLOT^RECIBO_OK_T.PRO\\(ONELKISTOGRAM)\\l  e-Nov- 1098  15:54 


«0-9T  9661-AON-9T\\(KVa00iSHr3N0)V\0adl“H0“0S133aV~X01dVpo3“oqp3jeVAsdpa\:|\\jBdonod\VKa 


E 


u(q  u|  uon«indod 
JO  uonovjj 


U|q  u]  uops(ndod 
JC  uonovjj 


?2  5 


-  I 
-  ' 
o|  £ 

oS  c 
o'* 


aiq  ui  aon«indod 

JO  uono«4i 


aiq  U]  aon«indod 
JO  uonsvJJ 


a 

“5  c 


uiq  u]  aoncindod 

JO  uonovJJ 


«*  .r 


2? 
I  E 


JO 


u}q  U(  uopvindod 
JO  uonooJtJ 


£ 

Si 

C 


<0  E 


2  g 


uiq  u|  uonvindod 

JO  uona«-iJ 


u|q  u]  uonvindod 
JO  uonooJJ 


u]q  u]  uon^indod 
JO  uonas^fj 


il  CL 

3  o  & 
I  ?  c 


JO  uonottjj 


ufq  U|  uon«indod 
JO  uonoojj 


o' 


o 


U{q  u]  uoneindod 
JO  uonoQJJ 


uiq  UI  uonutndod 
JO  uonooJJ 


uiq  UI  uopuindod 

JO  UOpOQJJ 


u|q  UI  uopaindod 
JO  uopavJJ 


uiq  UI  uopuindod 

JO  uonooJtj 


o  o  o  «o  o 

o  OI  M 

o  o  o  o  o 

o  d  d  o'  d  o'  o 

U|q  UI  uopuindod 
JO  UOpOQJtJ 


uiq  U|  uopoindod 
JO  uopoojj 


60 


Table  1.  Statistical  summary  of  linear-H  fits  to  Arecibo  World  Days  between  1985-1995. 
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Table  1.  (continues)  Statistical  sunimary  of  linear-H  fits  to  Arecibo  World  Days  between  1985-1995 
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Table  1.  (continues)  Statistical  summary  of  linear-H  fits  to  Arecibo  World  Days  between  1985-1995 
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Table  1.  Statistical  summary  of  linear-H  fits  to  Arecibo  World  Days  between  1985-1995 
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year,  B  =  sample  of  9801  fits,  C=sample  of  14000  profiles,  D=downward  maximum,  a  =  standard  deviation  of  sample. 


